Exact analytical results for quantum walks on star graph 
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In this paper, we study coherent exciton transport of continuous-time quantum walks on star 
graph. Exact analytical results of the transition probabilities are obtained by means of the Gram- 
Schmidt orthonormalization of the eigenstates. Our results show that the coherent exciton transport 
displays perfect revivals and strong localization on the initial node. When the initial excitation starts 
at the central node, the transport on star graph is equivalent to the transport on a complete graph 
of the same size. 
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CN| ■ The problem of coherent and non-coherent transport modeled by random walks has attracted much attention 
^ | in many distinct fields, ranging from polymer physics to biological physics, from solid state physics to quantum 
computation [l], Q, H, H|. Random walk is related to the diffusion models and a fundamental topic in discussions 
of Markov processes. Several properties of random walks, including dispersal distributions, first-passage times and 
encounter rates, have been extensively studied As a natural extension to the quantum world of the ubiquitous 

classical random walks, quantum walks (QWs) have also been introduced and widely investigated in the literature 0]. 
An important application of quantum walks is that QWs can be used to design highly efficient quantum algorithms. 
'■^ | For example, Grovers algorithm can be combined with quantum walks in a quantum algorithm for glued trees which 
provides even an exponential speed up over classical methods [8|, |9| . Besides their important applications in quantum 
algorithms, quantum walks are also used to model the coherent exciton transport in solid state physics [101 ]. It is 
■ shown that the dramatic nonclassical behavior of quantum walks can be attributed to quantum coherence, which does 
^ ' not exist in the classical random walks. 

There are two main types of quantum walks: continuous-time and discrete-time quantum walks. The main difference 
between them is that discrete-time walks require a coin which is just any unitary matrixplus an extra Hilbert space on 
7-H ■ which the coin acts, while continuous-time walks do not need this extra Hilbert space [111 ]. Apart from this difference, 
the two types of quantum walks are analogous to continuous-time and discrete-time random walks in the classical 
case (llj . Discrete-time quantum walks evolve by the application of a unitary evolution operator at discrete-time 
intervals, and continuous-time quantum walks evolve under a time-independent Hamiltonian. Unlike the classical 
case, the discrete-time and continuous-time quantum walks cannot be simply related to each other by taking a limit 
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as the time step goes to zero j!2j . Both the two types of quantum walks have been defined and studied on discrete 
structures in the market. 

Here, we focus on continuous-time quantum walks (CTQWs). Most of previous studies related to CTQWs in the 
. last decade concentrates on regular structures such as lattices, and most of the common wisdom concerning them 
relies on the results obtained in this particular geometry. Exact analytical results for CTQWs are also found for some 



particular regular structures, such as the cycle graph [13[ and Cayley tree [14[. Exact analytical results are difficult 



to get due to the cumbersome analytical investigations of the eigenvalues and eigenstates of the Hamiltonian. 

In this paper, we consider CTQWs on star graphs, and get exact analytical results for the first time. The star 
graph is one of the most regular structures in the graph theory and represents the local tree structure of the irregular 
and complex graphs. In mathematical language, a star graph of size N consists of one central node and N — 1 leaf 
nodes. All the leaf nodes connect to the central node, and there is no connection between the leaf nodes. Therefore, 
the central node has N — 1 bonds and the leaf nodes have only one bond. As we will shown, for such simple topology, 
we are able to derive exact analytical results for the transition probabilities. These analytical results exactly agree 
with the numerical results obtained by diagonalizing the Hamiltonian H using the software MATLAB. 

The coherent exciton transport on a connected network is modeled by the continuous-time quantum walks 
(CTQWs), which is obtained by replacing the Hamiltonian of the system by the classical transfer matrix, i.e., 
H = —T [lil Hjj |. The transfer matrix T relates to the Laplace matrix by T = —A. The Laplace matrix A has 
nondiagonal elements Aij equal to —1 if nodes i and j are connected and otherwise. The diagonal elements An 
equal to degree of node i, i.e., An — hi. The states \j) endowed with the node j of the network form a complete, 
ortho-normalised basis set, which span the whole accessible Hilbert space. The time evolution of a state \j) starting 
at time to is given by \j, t) = U(t, to)\j), where U(t, to) = exp[—iH(t — to)] is the quantum mechanical time evolution 
operator. The transition amplitude ctk,j(t) from state \j) at time to state \k) at time t reads ctk,j{t) — (k\U(t, 0)\j) 
and obeys Schrodingers equation [l7| . Then the classical and quantum transition probabilities to go from the state 
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\j) at time to the state \k) at time t are given by Pk,j(t) = (k\e~ tA \j) and 7r fcj (i) = \a k ^(t)\ 2 = \{k\e~ ttH \j)\ 2 0, 
respectively. Using E n and |q„) to represent the nth eigenvalue and ortho-normalized eigenvector of H, the classical 
and quantum transition probabilities between two nodes can be written as [l?], EH 

Pk,j(t) = J2 e ' tEn ( k M^\j), (i) 



^ 3 (t) = \a kj (t)\ 2 = \ 



(k\q„)(q n \j) 



(2) 



For finite networks, the classical transition probabilities approach the equal-partition l/N. However, the quantum 
transport does not lead to equal-partition. TTk,j(t) do not decay ad infinitum but at some time fluctuates about a 
constant value. This value is determined by the long time average of TTk,j(t) [l7l. fl8j. 



Xk,j 



= liniT^oo i J Q T 7Tfcj (t)dt 
= J2n.i( k \<ln)(qn\j)(j\qi){qi\k) 
xlimr-.ooi r n T e -«(^-^)^ 



(3) 
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where 5e„,Ei takes value 1 if E n equals to Ei and otherwise. In order to calculate Pk,j(t), ^kj(t) and Xk,j, all 
the eigenvalues E n and eigenstates |g n ) are required. For the star graph, in the following, we will first analytically 
calculate the eigenvalues and eigenstates, then give exact analytical results for the transition probabilities according 
to the above Equations. 

For a specific star graph of size N, we label the central node as node 1 while the leaf nodes are numbered as 2, 3, N. 
Because of centrosymmetric structure of the star graph, there are only four types of transition probabilities, namely, 
7Ti,i(i), 7T2,i(t) = 7Ti,2(i), Tt2,2(t) and 7T3,2(£) = 7T2,3(i). Transition probabilities between other nodes belong to these 
four types. Therefore, we only consider the four kinds of probabilities listed above. The Hamiltonian of the star graph 
can be written as, 
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H 



i=2 



The eigenvalues of the Hamiltonian have three discrete values: E\ = E% = ... = En -2 = 
En = N [19(. One set of eigenstates (i = 1,2, ...,N) corresponding to the eigenvalues is 



1, Eff-i 
\vi) = \i - 



(4) 

= and 

2) - |2) 



(i = 1,2, N — 2), \vff-i) = YliLi I*) an< ^ \ v n) = yii^i H) — N\l). However, this set of eigenstates are not orthogonal 
((^1^2) 7^ 0, etc), we use the Gram-Schmidt process [20| to orthogonalize this set of eigenstates. The Gram-Schmidt 
algorithm is a method for orthogonalizing a set of vectors in an inner product space 21| , the new orthogonal vectors 
— 1j 2, N — 2) are given by the following formula, 



where 



j=l v 3 ' J 1 

\v\) is applied in the iterative process 2^, 21|. According to the above equation, we obtain, 

2)-iEi- \i)< i= 1-2,..., N -2 



(5) 



l»- 

\Vi), 



j), i=l,2,...,N 
i = N-l,N. 



(0) 



The above new eigenstates {v'j}} is not normalized (i.e. 
the orthonormal basis {|<?i)} as follows, 




v^E-Lili), 



1, etc). After some algebraic calculations, we get 



(7) 
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One can easily prove is also a set of eigenstates of the Hamiltonian (i.e., H\qi) — Ei\qi)y i). {\qi}} forms 

an orthonormal and complete basis, satisfying (qi\qj) = &y and ^ - \qj)(qj\ = 1. Therefore, we can use this set of 
eigenstates {|<?i)} to calculate the transition probabilities in Eqs. |T]), ([2|) and ((3]). 

Substituting the orthonormal basis of Eq. into Eqs ([2]), we get iri,i(t) and 7T2,i(i) as follows, 

. . A 2 -27V + 2 2(JV-1) 



^.i(*) = J2- ^cosiVt. (9) 



For 7T2,2(i) and 7^2 (i), the calculation is analogous, but the expressions are cumbersome: 

7T2, 2 (t) = w ^ w [(N^-4N 3 + 5N 2 -2N + 2) 
+ (2N 3 - 6N 2 + AN) cos t 



2 **T\^,Kr ^ (!«) 



+ (2iV 2 -4iV) cos(A- l)t 
+(2A-2)cosAT*]. 



7T3,2(t) = nH n-i)'A (N 2 ~N+1) 

+ (N - N 2 )cost- Ncos{N~ l)t (11) 
+{N- l)cosJVi]. 

Other transition probabilities can also calculated, but they have the same expressions as Eqs (8)~(11). For instance, 
7T4 )2 (^) has the same analytical form as ^3.2, which is consistent with our intuition. Analogously, we can get the long 
time averages of the transition probabilities: 

X11 = (N 2 -2N + 2)/N 2 
X2,i = 2/iV 2 

X2,2 = (A^ 4 - 4iV 3 + 5N 2 -2N + 2) /N 2 /(N - l) 2 [ ' 

X3 ,2 = 2(N 2 - N + 1)/N 2 /{N - l) 2 . 

The transition probabilities depend on the size of the graph. In the thermodynamic limit of infinite network N — > 00, 
the transport displays high localizations on the initial position, i.e., Xi,j ~ O n the contrary, for the classical 

transport modeled by continuous-time random walks, the transition probabilities do not show any oscillation and 
approach to the equal-partition l/N at long times [19(. According to Eq. {1}, the four type transition probabilities 
can be written as, 

p 1(1 (t) = l/N + e- m (N-l)/N 
P2,i{t) = l/N-e- m /N 

p 2<2 (t) = l/N + (N — 2)/{N - l)e- 4 + e- Nt /N/(N - 1) [l6 > 
p 3 2 (t) = l/N - e-*/(N - 1) + e- Nt /N/{N - 1). 



In order to test the analytical predictions, we compare the classical Pk,j(t) predicted by Eq. (|13p with the numerical 
results obtained by numerically diagonalizing the Hamiltonian. The results for a star graph of N = 100 are shown in 
Fig. [U As we can see, the numerical results exactly agree with the analytical prediction in Eq. (fTB")) . The transport 
reaches the equal-partitioned distribution l/N at long times. However, when the excitation starts at the central node, 
the transport reaches the equal-partition more quickly than the transport which starts at the leaf nodes (Compare 
the curves in Fig. [T]). 

For the quantum transport, we compare the transition probabilities in Fig. [5J The numerical results (marked as 
black squares) exactly agree with the theoretical results in Eqs. (|8"j)- ([TTj) . We note that all the transition probabilities 
show periodic recurrences. Comparing iri,i(t) and 1^2,2^) (See Fig. [2] (a) and (c)), we find that there are high 
probabilities to find the exciton at the initial node. This suggests that the coherent transport shows high localizations 
on the initial nodes [l9| . The oscillation amplitudes of the return probabilities 7ri j i(i) and 7T2,2(i) are comparable but 
the oscillation periods are quite different. The oscillating period of ^2,2^) is 100 (N) times of that of 7ri,i(<). This 
could be interpreted by the analytical expressions in Eqs. (jHJl and (| 10[) . Similar behavior also holds for 7T2,i(t) and 
7T3, 2 (i) (See Fig. [5](b) and (d)), but the oscillation amplitude is smaller than the return probabilities. This also can be 
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FIG. 1: (Color online) Classical transition probabilities Pk,j(t) versus t. The marked points are numerical results and the curves 
are analytical predictions in Eq. (|13|l . 




FIG. 2: (Color online) Quantum transition probabilities 7ri.i(i) (a), ■K2,i{t) (b), 7T2,2(i) (c) and 7r3,2(£) (d). The points denoted 
by black squares are numerical results and the curves are theoretical predictions in Eqs. (|8])- (|lip . 

understood from the analytical results in Eqs. ^ and (fTTj) , where the transition probability is mainly determined by 
the high order term of N. The small value of oscillating period of 7r 2i i(t) and Hi ,i(i) suggests that there are frequent 
revivals when the exciton starts at the central node, compared to the transport starting at the leaf nodes. 

The quantum limiting probabilities in Eq. (|12p are only a function of graph size N . Fig. [3] shows the quantum 
limiting probabilities for numerical results and theoretical predictions. Both the results agree with each other. We find 
that the return probabilities and %2,2 are an incremental function of TV and approach to 1 in the limit N — ► oo. 
By contraries, %2,i and X3,2 decrease with N and close to 2/N 2 in the limit N — ► oo. We note that Xl,i differs from 
X2,2 for small values of N. Such deviation diminishes as N increases. This suggests that the strength of localizations 
is almost the same for central-node and leaf-node excitations. The only difference is that the frequency of revivals 
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FIG. 3: (Color online) Long-time limiting probabilities XMj X2,2 (a) and X2,i, X3,2 (b) as a function of the network size N. The 
points marked as symbols are numerical results and the curves are analytical results predicted by Eq. (|12|) . 



(oscillation period) for central node excitation is much higher than that for leaf node excitation. 

To address the similarity and difference between the star graph and complete graph, we proceed to consider the 
transport on a complete graph of size N. The complete graph is fully connected, thus the Hamiltonian is given by 
H = (N — l)J2i — Si^j KKil- The eigenvalues are two different values: E\ = E 2 = ... = £jv-i = N and 
En — 0. One set of un-orthogonal states corresponding to the eigenvalues can be written as: \vi) = \i + 1) — |1) 

(i = 1, 2, N— 1) and \vn) = zCj=i li)- Using the Gram-Schmidt orthonormalization (See Eq. the orthonormal 
basis for a complete graph is, 




Substituting the above Equation into Eq. ([2]), we get the quantum transition probabilities for complete graph, 

f jV 2 -27V+2 , 2(N-1) m . _ ■ 

i^j. (15) 

We note that Eq. (|T5|) is exactly the same form as Eqs. ([5]) and ©. This indicates that the transport starting at 
central node on star graph is equivalent to the transport on a complete graph of the same size. 

In summary, we have studied coherent exciton transport of continuous-time quantum walks on star graph. Exact 
analytical results of the transition probabilities are obtained in terms of the Gram-Schmidt orthonormalization. We 
find that the coherent transport shows perfect recurrences and there are high frequency of revivals for central node 
excitation. Study of long time averages suggests that the quantum transport displays strong localizations on the 
initial node. When the initial excitation starts at the central node, the transport on star graph is equivalent to the 
transport on a complete graph of the same size. 

This work is supported by National Natural Science Foundation of China under projects 10575042, 10775058 and 
MOE of China under contract number IRT0624 (CCNU). 
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